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ABSTRACT 

I I present an analysis of the XMM-Newton observations of four millisecond 

§ : pulsars, J0437-4715, J2124-3358, J1024-0719, and J0034-0534. The new data 

^ ' provide strong evidence of thermal emission in the X-ray flux detected from the 

first three objects. This thermal component is best interpreted as radiation from 
pulsar polar caps covered with a nonmagnetic hydrogen atmosphere. A nonther- 
mal power-law component, dominating at energies E >3 keV, can also be present 
in the detected X-ray emission. For PSR J0437-4715, the timing analysis reveals 
> ; that the shape and pulsed fraction of the pulsar light curves are energy dependent. 

This, together with the results obtained from the phase-resolved spectroscopy, 
CN ' supports the two-component (thermal plus nonthermal) interpretation of the pul- 

Q ■ sar's X-ray radiation. Highly significant pulsations have been found in the X-ray 
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flux of PSRs J2124-3358 and J1024-0719. For PSR J0034-0534, a possible X-ray 
counterpart of the radio pulsar has been suggested. The inferred properties of the 



CL|' detected thermal emission are compared with predictions of radio pulsar models. 
O 

^ ; Subject headings: pulsars: individual (J0437-4715, J2124-3358, J1024-0719, 

^ ■ J0034-0534) — stars: neutron — X-rays: stars 
> 

! 1. Introduction 

Of more than 1500 currently known rotation-powered (radio) pulsars^, 122 possess very 
short spin periods P < 10 ms. These millisecond pulsars are generally believed to be very 
old neutron stars (NSs), spun up by accretion in binary systems, with characteristic ages 
r = P/2P ~ 0.1-10 Gyr and low surface magnetic fields Esurf oc (PP)^/^ _ iqS-IO^o G 
(e.g., Taylor et al. 1993). As millisecond pulsars are intrinsically faint at optical wavelengths 
and most of them reside in binary systems with optically brighter white dwarf companions. 



^Space Science Laboratory, NASA MSFC SD50, Huntsville, AL 35805; vyacheslav.zavlin@msfc.nasa.gov 

^According to the catalog provided by the Australia Telescope National Facility (ATNF), 
http://www.atnf.csiro.au/research/pulsar, for the middle of 2005. 
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X-ray energy range is the main source of information on these objects outside the radio band. 
So far firm X-ray detections have been reported for nine miUisecond (non-variable) pulsars, 
either solitary or in binaries without interaction with companions. They are given in Table 1 
(this list does not include X-ray emitting millisecond pulsars in the global clusters 47 Tuc 
and NGC 6397 because of a lack of a detailed information on their X-ray properties — see, 
e.g., Grindlay et al. 2002). 

X-ray emission from radio pulsars consists of two different components, thermal and 
nonthermal (see, e.g., Pavlov et al. 2002, and Kaspi et al. 2006 for reviews). The nonthermal 
component is described by a power-law (PL) spectral model and attributed to radiation 
emitted from pulsar magnetosphere, whereas the thermal emission can originate from cither 
the whole surface of a coohng NS or small hot spots around the magnetic poles (polar caps; 
PCs) on the star surface, or both. As predicted by virtually all pulsar models, these PCs 
can be heated up to X-ray temperatures (~ 1 MK) by relativistic particles generated in 
pulsar acceleration zones. In the case of millisecond pulsars, the entire surface at a NS 
age of ~ 1 Gyr is too cold, <0.1 MK, to be detectable in X-rays (although it may be 
seen in UV/FUV — see §2). Therefore, only nonthermal and/or thermal PC radiation is 
expected to be observed in X-rays from these objects. The pulsars presented in Table 1 
can be divided in two distinct groups (see also Kuiper & Hermsen 2004). The first group 
consists of PSRs B1937+21, B1957+20, J0218+4232, and B1821-24, with high rotational 
energy losses, E > 10^^ erg s~^, and estimated X-ray luminosities, Lx > 10^^ erg s~^. Their 
X-ray emission is of a nonthermal origin, with pure PL spectra of photon index P = 1.5- 
2 and pulsed profiles with strong narrow peaks and large pulsed fractions, /p > 50% (see 
references in Table 1). Probably, PSR J0751-I-1807 can also be associated with this group 
because its spectrum is better fitted with a PL model of P ^ 1.6. Three other pulsars, 
J0437-4715, J0030+0451 and J2124-3358, with lower spin-down energies, E < 10^^ erg s-\ 
belong to the second group. X-ray emission from these NSs shows significant contribution of 
a thermal component, with smoother pulsations and lower pulsed fractions. The properties of 
X-rays detected from PSR J1012-I-5307 are rather uncertain. However, there is an indication 
of a thermal component in its X-ray flux. Current theoretical models are not elaborate 
enough to predict in which millisecond pulsars the thermal PC component would prevail 
over the nonthermal one. It, however, seems certain that both the thermal and nonthermal 
luminosities should increase with rotational energy loss E. Therefore, only analysis of X-ray 
observations can help in elucidating the radiative properties of millisecond pulsars. 

The high sensitivity and spectral resolution of Chandra and XMM-Newton provide new 
opportunities for studying X-ray emission from faint celestial objects. This paper presents an 
analysis and interpretation of new XMM-Newton observations of four millisecond pulsars, 
J0437-4715, J2124-3358, J1024-0719, and J0034-0534 (§§2-5). The obtained results are 
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summarized and discussed in § 6. 

2. PSR J0437-4715 

This is the nearest and X-ray brightest miUisecond pulsar known. It is in a 5.7 d binary 
orbit with a low-mass white dwarf companion. Its pulsed X-ray emission was first detected 
with ROSAT in the 0.1-2 keV range (Becker & Triimper 1993). Those data could be fitted 
with either a single PL model of F ^ 2.5, suggesting a nonthermal origin of the detected 
emission, or a thermal PC model consisting of two opposite PCs covered with nonmagnetic 
hydrogen atmosphere'^ and temperature decreasing from ~ 1.5 MK in the PC "core" of 
about ~ 0.3 km radius to ~ 0.4 MK in the PC "rim" of a ~ 3 km radius (Zavlin & Pavlov 
1998). Becker & Triimper (1993) found that pulsed fraction /p of the pulsar flux increases 
with photon energy E from about 30% at ~ 0.2 keV up to about 50% at 1 keV (although 
Becker & Triimper 1999 found no energy dependence of /p from a later ROSAT observation 
of PSR J0437-4715). On the other hand, thermal emission from NS atmospheres is essentially 
anisotropic, even at low Bsuvi, with radiation intensities beamed stronger at higher E (Zavlin 
et al. 1996), that can result in an increase of fp with photon energy. Therefore, if the energy 
dependence of fp found by Becker & Triimper (1993) is real, it could be regarded as an 
additional indication of the thermal origin of the detected radiation. The pulsar was also 
observed with Chandra. Analysis of the combined ROSAT and Chandra data (Zavlin et 
al. 2002) showed that the 0.1-7 keV spectrum of PSR J0437-4715 cannot be described with 
a simple one-component model and requires two components: a nonthermal PL spectrum 
of r f=i 2.0 and a thermal PC model similar to that suggested by Zavlin & Pavlov (1998) 
from the ROSAT data alone. The thermal (bolometric) and nonthermal (in 0.2-10 keV) 
luminosities were found to be Lboi ~ 2.3 x lO^'^ erg s~^ and f« 0.6 x 10^° erg s~^ (for 

a distance d = 0.14 kpc — see Table 1), respectively. In this interpretation, the nonthermal 
component dominates at energies E >3 keV. The Chandra, data provided also an energy- 
integrated pulse profile with /p ^ 40% and the X-ray peak at the same phase as the radio 
peak. In addition, PSR J0437-4715 is the only millisecond pulsar detected near the optical, 
at the wavelengths A = 1155-1702 (Kargaltsev et al. 2004), that provided an estimate on 
the temperature of the entire NS surface, Tgurf ~ 0.1 MK — surprisingly high for such an 
old NS, implying special heating mechanism(s) to operate in this pulsar. 



^See Zavlin & Pavlov (2002) for a recent review on NS atmosphere models. 
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2.1. XMM-Newton observation 

PSR J0437-4715 was observed with XMM-Newton on 2002 October 9 (orbit 519) for 68.3 
and 67.2 ks effective exposures with the EPIC^-MOS and EPIC-pn instruments, respectively. 
EPIC-MOSl and M0S2 were operated in Full Frame mode providing an image of an area of 
a ~ 14' radius and a time resolution of 2.6 s. EPIC-pn was in Timing mode. In this mode 
only one EPIC-pn CCD is in use and the EPIC-pn image is collapsed in one-dimensional 
(1-D) count distribution in the direction perpendicular to the CCD read-out, that provides a 
0.03 ms time resolution suitable for detecting pulsations in X-ray flux of millisecond pulsars. 
Thin filters were used for each EPIC instrument. The EPIC data were processed with the 
most recent SAS package^ (ver. 6.1.0). 

Figure 1 presents a 20' x 20' image constructed from combined EPIC-MOSl and M0S2 

data. PSR J0437-4715 is the second brightest source in the filed of view. The EPIC-pn 
field of view is also shown in Figure 1, as well as the direction of the CCD read-out. It is 
evident that the satellite position angle (138.7°) chosen for observing PSR J0437-4715 was 
not optimal, which led to a partial contamination of the pulsar data in the 1-D EPIC-pn 
image by the much brighter neighboring AGN. This resulted in an asymmetric shape of the 
1-D distribution of EPIC-pn counts shown in the upper panel of Figure 2. To evaluate the 
contamination, we constructed three 1-D images from the EPIC-MOS data confined within 
the EPIC-pn field of view with the same column width of 4'.'2 as in the EPIC-pn CCD: 
from the areas above and below the short-dashed line in Figure 1, which contain negligible 
fractions of pulsar and AGN quanta (respectively), and from the total area. Examination of 
these three distributions plotted in the lower panel of Figure 2 shows that the contamination 
of the pulsar is insignificant in the CCD columns below # 40, which were used for the analysis 
below. 



2.2. Spectral analysis 

Both EPIC-MOS spectra of PSR J0437-4715 were extracted from circles of a 30" radius 
centered at the pulsar position, each containing about 82% of all detected source counts. 
Background was estimated from a few source-free regions. The extracted spectra were binned 
with at least 30 source counts per bin. The estimated total count rates of the pulsar are 
161 ± 2 and 167 ± 2 counts ks"^ in the 0.3-10 keV range for the EPIC-MOSl and M0S2 



^European Photon Imaging Camera 
^http : //xmm . vilspa . esa . es 
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detectors, respectively (energies below 0.3 keV were not used because of uncertainties in 
the EPIC-MOS calibration). The EPIC-pn spectrum was extracted from the CCD columns 
#28-39, what provides about 85% of all detected source counts. Background was estimated 
from the columns 7^ 12-23. The EPIC-pn spectrum was binned with at least 100 source 
counts in each bin. The total source count rate was measured to be 674 ±5 counts ks~^ in the 
0.3-10 keV range (there is no spectral information below 0.3 keV in the EPIC-pn detector in 
Timing mode) . All instrument response matrices and effective areas were generated with the 
SAS tools rmfgen and arfgen (these tools were also used in analysis of all XMM-Newton data 
discussed in this paper). The source count rates in the Reflection Grating Spectrometers, 
~ 20 counts ks~^, are too low for a meaningful analysis of this data. 

Applying a simple one-component model to the three combined EPIC spectra resulted 
in unreasonable model parameters and unacceptable flt quality. A best single (absorbed) 
PL flt results in a photon index of F ^ 3.5 and absorbing hydrogen column density of 
?T'H,2o — nH/(10^'^cm^) ^ 11.6 (vs. the estimate nH,2o = 0.1-0.3 obtained from independent 
measurements — see Zavlin & Pavlov 1998), with the minimum value of xl = 2.47 (for 
1/ = 426 degrees of freedom). Even worse fit is yielded by a one-temperature thermal model, 
6-g-j xl = 7.14 for a single blackbody (BB) spectrum, because of a large data excess at 
energies E >2 keV. A two-component model, BB-I-PL, results in xl = 1-67 (t/ = 424). With 
account for problems in the cross-cahbration of the EPIC instruments^ it can be considered 
as a reasonably good flt. In this model the PL component with P = 2.9 provides about 
70% of the observed X-ray flux {fx ~ 1.1 x 10~^^ erg s~^ cm~^ in 0.2-10 keV). However, 
the obtained value of nH,2o = 1-4 still signiflcantly exceeds the independent estimates on the 
hydrogen column density. In addition, the inferred photon index is much larger than P = 
1.4 — 2 found in nonthermal emission from other radio (including millisecond ones) pulsars. 
A similar problem arises when one-temperature PC model with hydrogen atmosphere is 
applied. Hence, we regard this model as rather infeasible. The next option is a PC model 
with a nonuniform temperature distribution combined with a PL component. Applying the 
"core" -I- "rim" model for two PCs covered with a hydrogen atmosphere (Zavhn et al. 1996) 
and assuming standard NS mass M = 1.4Mq and radius = 10 km yields''' T^°^^ = 1.4 ±0.2 
MK and Tp™ = 0.52 ± 0.16 MK (unredshifted values), = 0.35 ± 0.18 km and i?™ = 
2.6±0.4 km, and P = 2.0±0.4, for nH,2o varying between 0.1 and 0.3, with the minimal value 
of xl — 1-57 {u — 422). Formally, according to an F-test, the probabihty that the spectral 



^See Kargaltscv et al. (2005) for results of the spectral analysis of XMM- 
Newton data on the NS RX J1856. 5-3754, as well as other examples at 
http : //xmm . vilspa . esa . es/docs/ docuiiients/CAL-TN-0018-2-4 . pdf . 

''Errors in fitting parameters are given at a Icr level for one interesting parameter. 
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data on PSR J0437-4715 require a two-temperature PC model (instead of a one-temperature 
thermal component) is 99.99992%. This model implies equal angles ( (between the line-of- 
sight and pulsar spin axis) and a (between the magnetic and spin axes) of 45°. Figure 3 
presents a best fit with this two-component model. The estimated bolometric luminosity 
of two PCs is Lboi = (3.4 ± 0.7) x 10^° erg s"\ and the PL luminosity in 0.2-10 keV is 
2^nonth _ (^0.5 ±0.2) X 10^° crg s~^. Changing the pulsar geometry (angles ( and a) as well as 
the NS mass-to-radius ratio in a plausible range M/R = (1.0-1.8) M0/(1O km) responsible 
for the gravitational effects (redshift and bending of photon trajectories) increases the errors 
in the PC parameters, but insignificantly it affects the inferred thermal luminosity (see Zavlin 
& Pavlov 2004 for more details). Generally, the spectral results obtained from the EPIC 
data are in good agreement with those derived from the previous observations. 

Another hypothesis is that the pulsar spectrum is of a pure nonthermal origin with the 
PL slope changing somewhere in the X-ray range. To verify this, a broken PL model was 
applied to the EPIC spectra, that yielded a best fit with nn 20 ~ 2.3±0.5, and photon indices 
Pi = 2.4 ± 0.2 and T2 = 3.6 ± 0.1 below and above the break energy £^br = 1-05 ± 0.05 keV 
{xt = 1-73 for u = 424). However, this model significantly overpredicts the UV/FUV fiuxes 
detected from PSR J0437-4715 (Kargaltsev et al. 2004). This, together with the large value 
of the interstellar absorption inferred in the fit, makes the broken-PL interpretation rather 
implausible. 



2.3. Timing analysis 

For the timing analysis, we used 38393 counts extracted from the EPIC-pn CCD columns 
7^35-39 in the 0.3-6 kcV range chosen to maximize the signal-to-noise ratio {S/N = 163). 
Of those counts, 88% belong to the pulsar. The photon arrival times were transformed 
to the solar system barycenter with the barycen tool of the SAS package (this procedure 
was also applied to all other timing data discussed in this paper) . As timing parameters of 
millisecond pulsars are known to be very stable, we invoked the binary ephemeris parameters 
from Zavlin et al. (2002). These parameters used in the test {n is number of harmonics 
— see Buccheri et al. 1983) yielded the most significant value = 2170.2 with n = 2. The 
maximum of .Z'ljj^^^^ = 2199.1 is reached at a frequency / which differs from the radio value 
by Sf = 0.8 yuHz, corresponding to a relative accuracy, Sf/f ~ 5 x 10~^, consistent with 
that estimated for the EPIC-pn timing^. 

We extracted the pulsed profiles of PSR J0437-4715 in five energy ranges, 0.3-0.5, 0.5- 



See http : //xmin . vilspa . esa . es/docs/documents 
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0.8, 0.8-2, 0.3-2 keV, and 2-6 keV, which are given in Figure 4 together with estimated values 
of pulsed fraction /p. The statistical test described in Zavlin & Pavlov (1998) shows that 
the shapes of the three lower-energy light curves are actually the same (the probability that 
they are different is less than 90%). However, the 2-6 kcV light curve seems to be more 
symmetric, with a peak apparently narrower than those in the lower-energy pulse profiles. 
One may also speculate that there is a phase shift of 50 ~ 0.1 between the peaks at energies 
below and above 2 keV. This apparent difference is supported by the same statistical test 
yielding a fairly high probability of 99.99997% (or a 4.2(7 significance) that the shapes of 
the 0.3-2 and 2-6 keV pulse profiles are different. Another feature in the extracted fight 
curves is the change of the source intrinsic pulsed fraction with energy: /p increases by 
about 10% from ~ 32% at the lowest energies to ~ 42% at ~ 1 keV, similar to the 
energy behavior of /p first found by Becker & Triimper (1993) from the ROSAT data. The 
estimated pulsed fraction at £' > 2 keV is less certain because of a scantier statistics and 
strong background contamination at higher energies, although it may be as high as ~ 65%. 
Generally, the energy dependence of the pulsar light curves can be attributed to different 
emission mechanisms producing the X-ray fiux of PSR J0437-4715, as suggested by the 
thermal-plus-nonthermal model. 

Figure 5 shows the pulse profiles of PSR J0437-4715 obtained from the ROSAT Position 
Sensitive Proportional Counter, Chandra High Resolution Camera (see Zavfin et al. 2002) 
and XMM-Newton EPIC-pn data (in 0.3-6 keV). The pulsed fractions in the latter two 
fight curves are similar to each other and somewhat larger than that found in the ROSAT 
data. This may be explained by both the properties of the pulsar emission (increase of /p 
with energy) and the greater sensitivity of the XMM-Newton and Chandra instruments to 
higher-energy photons. The shapes of these three light curves are clearly asymmetric, with a 
longer rise. This asymmetry could be caused by contribution of the nonthermal component 
(more than 10% at all energies) whose peak is shifted in phase with respect to the pulse of 
the thermal emission. 

As the error in the EPlC-pn absolute timing can be as large as 0.5 ms, the phasing 
between the radio and X-ray peaks was not attempted (this is also one of the reasons why 
the zero phases in Fig. 5 are arbitrary). 

2.4. Phase-resolved spectroscopy 

In principle, the EPIC-pn data on PSR J0437-4715 allow one to perform a phase- 
resolved spectroscopy of the pulsar emission. To do this, we extracted the pulsar spectra 
in five equal intervals between the rotational phases and 1, with the zero- phase as in 
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Figures 4 and 5, from the same CCD columns chosen for the phase-integrated spectrum 
(§ 2.2). The spectra were binned to collect at least 50 source counts per bin. A simplest and 
broadly used approach is to fit the phase-resolved spectra with the same model to determine 
the dependences of model parameters on rotational phase. Based on the results of §2.2, 
we applied the two-temperature PC ( "core" + "rim" ) plus PL model with temperatures of 
the thermal components and photon index fixed at the best values derived for the phase- 
integrated spectrum. The hydrogen column density was also fixed at nH,2o ~ 0.2. Figure 6 
shows phase dependences of the three fitting parameters, sizes of the PCs and nonthermal 
fiux, normalized by the corresponding values found in the phase-integrated spectral analysis. 
The former two represent the PC areas projected onto a plane perpendicular to the line-of- 
sight of a distant observer. As it could be expected, these two dependences are very similar 
to each other and resemble the shape of the lower-energy light curves {dX E <2 keV) shown 
in Figure 4. On the other hand, the phase dependence of the nonthermal flux clearly differs 
from those of the PC areas, with a narrower peak shifted with respect to the maxima in the 
PC phase dependences, indicating that the shift between the peaks in the pulsed proflles 
of PSR J0437-4715 extracted at energies below and above 2 keV (Fig. 4) is real. Also, the 
modulation of the nonthermal fiux is rather large, by about 60%, pointing to a high intrinsic 
pulsed fraction of the nonthermal component. 

More detailed information could be obtained from combining the phase-resolved spec- 
troscopy with modeling the pulsar light curves. However, such an analysis is beyond the 
scope of this paper. 

3. PSR J2124-3358 

X-ray emission from this solitary pulsar was first detected with ROSAT (Becker & 
Triimper 1999). These data provided a 4(T detection of pulsations of the pulsar X-ray fiux 
with a pulsed fraction /p ~ 33%. PSR J2124-3358 was also observed with ASCA (Sakurai 
et al. 2001). No significant pulsations were found in those data of a scanty statistics. The 
pulsar spectrum could be equally well fitted with either a single PL spectrum of photon 
index P ~ 2.8 or a thermal BB model of a temperature^ ~ 3.6 MK emitted from an area 
of an apparent radius ~ 0.02 km (for a distance d = 0.27 kpc). Based on the inferred 
model parameters, Sakurai et al. (2001) concluded that the thermal interpretation is more 
plausible. 



^The superscript "cx)" stays for redshifted values as measured by a distant observer. 
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3.1. XMM-Newton observation 

XMM-Newton observed PSR J2124-3358 on 2002 April 14-15 (orbit 430) for 68.9 and 
66.8 ks effective exposure for the EPIC-MOS and EPIC-pn detectors operated with medium 
and thin filters, respectively, in the same observational modes as for PSR J0437-4715. Fig- 
ure 7 shows a 20' x 20' combined EPIC-MOSl and M0S2 image with the pulsar located close 
to the image center. With the sateUite position angle of 76.7° in this observation, the pulsar 
emission detected in the EPIC-pn data was contaminated by weaker field sources (see Fig. 7 
for the EPIC-pn field of view) . Analysis of the 1-D distribution of EPIC-pn photons shows 
that the contamination is substantial only at energies E >2 keV. Hence, the 0.3-2 keV range 
was used for the analysis of the pulsar radiation detected with the EPIC-pn instrument. 



3.2. Spectral analysis 

Two EPIC-MOS spectra PSR J2124-3358 were extracted from circles of a 40" radius 
centered at the pulsar position, encircling about 88% of all detected source counts. The 
extracted spectra were binned with at least 25 source counts per bin. It was found that the 
pulsar emission is completely buried under background at energies above 3 keV. Hence, only 
the 0.3-3 keV range was used for the EPIC-MOS data on the pulsar (0.3 keV is the calibration 
"threshold" for these instruments — see § 2.2). The inferred total source count rates are 13±2 
and 17 ± 2 counts ks"^ for the EPIC-MOSl and M0S2 detectors, respectively. The EPIC- 
pn source (plus background) spectrum was extracted from the CCD columns # 35-42 (see 
insert in Fig. 7) . This was estimated to provide about 82% of all pulsar counts detected with 
EPIC-pn. Background was evaluated from the columns # 51-58. The EPIC-pn spectrum of 
PSR J2 124-3358 in the 0.3-2 keV range was binned with at least 40 source counts in each 
bin. The total source count rate in this instrument is 57 ± 3 counts ks~^. 

Estimates for the hydrogen column density in the direction to PSR J2124-3358 are 
rather uncertain. The pulsar dispersion measure DM = 4.6 pc cm~^ pc suggests nH,2o ~ 1, 
whereas estimates^° for objects close to PSR J2124-3358 indicate on larger values, nH,20 ~ 3. 
The total Galactic column density in the direction to PSR J2124-3358 is nni — 6 x 10^° 
cm~^. Hence, nH,2o = 1^3 can be regarded as a plausible range. The three EPIC spectra 
were first fitted with a single PL model. Although this fit is formally acceptable {xl = 1-1 for 
1/ — 121), it yielded a very large photon index P = 3.3 ± 0.4 and too high hydrogen column 
density nH,2o = 16 ± 5. These two parameters make the pure nonthermal interpretation 



See http: //archive . stsci . edu/euve/ism/ismf orm.html 
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of the pulsar X-ray emission infeasible. Applying a single thermal model, either BB or 
onc-tcmpcrature PCs with hydrogen atmosphere leaves a significant data excess at energies 
E >1.5 keV, suggesting that one more component is required to fit the spectra. 

First two-component model to test was one-temperature PCs with hydrogen atmosphere 
plus a PL spectrum. Such a model yielded a good fit with xl — 1-1 — 119) and hydrogen 
column density t7.h.2o = 3 ± 2 compatible with the independent estimates. As the pulsar 
emission is detected only below 3 keV, the slope of the nonthermal component is rather 
unconstrained, F = 2.1 ± 0.7. Nevertheless, the (isotropic) nonthermal luminosity is well 
determined, L'^°°*'^ = (0.9 ± 0.2) x 10^° erg s'^ in the 0.2-10 kcV range (for d = 0.27 kpc). 
The inferred PC parameters are Tpc = 1.3 ±0.1 MK and Rpc = 0.32 ± 0.04 km (for the 
standard NS mass M — IAMq and radius R — 10 km, and pulsar angles ( = a — 45°). 
The corresponding bolometric luminosity of two PCs, L^oi = (1-0 ± 0.2) x lO^'^ erg s~^, is 
almost the same as the estimated PL luminosity. A best PC-plus-PL fit is shown in Figure 8. 
As hydrogen atmosphere spectra are much harder than the BB ones at the same effective 
temperature, using the BB model for the thermal component results in a higher temperature, 

~ 2.4 MK, and a much smaller radius of emitting area, ~ 0.04 km (see, e.g., Zavlin 
& Pavlov 2004 for more examples). 

Next, a two-temperature PC ( "core" -|- "rim" ) model was applied to the EPIC spectra 
of PSR J2124-3358. It yielded a fit of the same quality as the one-temperature PC plus PL 
fit and the following model parameters: T^^'"" = 2.2 ± 0.2 MK, T™^^ = 0.5 ± 0.1 MK, and 
R'^J'^ = (0.11 ± 0.03) km, -R™ = (1.9 ± 0.7) km (for the same assumption on the NS mass, 
radius and geometry as above). Interestingly, these numbers are similar to those estimated 
for PSR J0437-4715 (§2.2). The bolometric luminosity of these PCs is Lboi = (l.SiO.l) xlO^^ 
erg s~^. No other spectral component (e.g., PL) in addition to this PC model is required 
to fit the observed data. An upper hmit on nonthermal emission of PSR J2 124-3358 is 
^^nonth ^ Q Qg ^ ;^q30 gj.g g-1 ^^le 0.2-10 keV range). 



3.3. Timing analysis 

To obtain a maximal signal-to-noise ratio {S/N = 21), the EPIC-pn CCD columns 
# 37-40 were used for extracting source (plus background) counts in the 0.3-2 keV range. 
Of total 8506 counts extracted, 39% was estimated to be emitted from the pulsar. Using 
the radio ephemeris parameters of PSR J2124-3358 from the ATNF catalog, spin frequency 
/o = 202.793897234988 Hz and its derivative / = -8.447 x 10"^^ s'^ (at MJD 50288.0), in 
the Zjj^ test immediately revealed a signal, with the most significant component Z| = 97.1. 
The probability to obtain such a value by chance in one trial is 4.1 x 10~^°. The maximum 
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value ^Imax — is reached at a frequency which differs from /o by 5f — 0.7 /xHz, similar 
to the case of PSR J0437-4715 (§2.3). 

Because of the strong background contamination and relatively narrow energy range 
where the pulsar data are available, only the energy- integrated (in 0.3-2 keV) light curve 
of PSR J2 124-3358 was extracted (see Fig. 9). It reveals one broad prominent pulse per 
period (at phase ^ 0.8 in Fig. 9) and a possible weaker peak separated from the main 
one by 6(f) ~ 0.35 (or 6(f) ~ 0.65, depending on which pulse is leading). Alternatively, the 
obtained pulse may be described as a single broad peak with an asymmetric shape, a steeper 
rise and a longer trail, as predicted by relativistic effects (in particular, the Doppler boost) 
in fast rotating pulsars (Braje et al. 2000). The estimated source intrinsic pulsed fraction, 
/p = 56 ± 14% is fairly high, but it is still consistent with thermal PC models, depending on 
the NS geometry and compactness (the M/R ratio). 

4. PSR J1024-0719 

This is a solitary pulsar with a spin period P ^ 5.2 ms, characteristic age r ~ 4.4 
Gyr and rotational energy loss E ~ 5.3 x 10^'^ erg s^^ (according to the ATNF catalog). It 
is a relatively close object located at a distance d ~ 0.39 kpc, as derived from the pulsar 
dispersion measure DM — 6.49 cm"^ pc and the model Galactic distribution of free electrons 
by Cordes & Lazio (2003). From ROSAT data Becker & Triimper (1999) suggested an X-ray 
counterpart for PSR J1024-0719 (at a 4a significance level). Because of very poor statistics 
available (~ 25 counts) and properties of the instrument used in the ROSAT observation, 
neither spectral nor timing information on the pulsar X-ray counterpart could be obtained. 

4.1. XMM-Newton observation 

PSR J1024-0719 was observed with XMM-Newton on 2003 December 2 (orbit 729) for 
68.0 and 66.2 ks effective exposures for the EPIC-MOS and EPIC-pn instruments, respec- 
tively, with the same filter and observational modes as for PSR J0437-4715. Figure 10 shows 
a 20' X 20' combined EPIC-MOSl and M0S2 image in the 0.3-2 keV range. PSR J1024-0719 
is clearly detected at a position which differs from its radio position by only 0'.'4, that is well 
within the 2"-3" uncertainty of the EPIC absolute astrometry. As seen in the EPIC-MOS 
image, the EPIC-pn field of view contains many background sources, with a comparable (or 
even larger) brightnesses, which heavily contaminate the 1-D image of PSR J1024-0719 in 
the EPIC-pn data (see insert in Fig. 10). This makes the EPIC-pn data virtually useless 
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for spectral analysis of the pulsar emission because no reliable background subtraction is 
possible. 

4.2. Spectral analysis 

Only EPIC-MOS data were used to evaluate spectral properties of PSR J1024-0719. 
The pulsar counts were extracted from circles of a 20" radius, which contain about 70% of 
all photons detected from PSR J1024-0719. Background was estimated from a few source- 
free regions. The extracted spectra were binned with at least 20 source counts per bin in the 
0.3-2 keV range and 9-11 source counts in the bins at 2-3 keV. The total source count rates 
are 2.3±0.3 and 2.7±0.3 counts ks"^ for the EPIC-MOSl and M0S2 detectors, respectively. 

The pulsar dispersion measure gives a hydrogen column density nH,2o ~ 2 toward 
PSR J1024-0719. Independent estimates of interstellar absorption for objects close to the 
pulsar suggest nH,20 ~ 6, whereas the total Galactic column density is nni — 5 x 10^° cm~^. A 
single PL fit to the pulsar spectra produced a best result with F = 3.7±0.8 and nH,20 = 23 ±5 
{x1 = 1.4 for u = 11). Similar to the case of PSR J2124-3358, these inferred parameters 
make the pure nonthermal interpretation hardly plausible. A single BB model fits the spec- 
tra equally well {xl = l-l)- The obtained parameters are nH,2o = 2 ± 2 (compatible with the 
independent estimates), = 2.9 ± 0.3 MK and = 0.03 ± 0.01 km (for d = 0.39 kpc). 
One-temperature PC model with nonmagnetic hydrogen atmosphere (and the same assump- 
tion on the NS parameters as in §§ 2.2 and 3.2) yielded Tpc = 1.8±0.4 MK and Rpc = O.liO.l 
km (xl = 1-0)) the same nH,2o range as in the BB fit. A best PC model is shown in 
Figure 11. The estimated bolometric luminosity of two PCs is Lboi = (0.4 ± 0.2) x lO^'' erg 
s~^. Adding a PL component of F = 2 to the best PC model fit put virtually the same upper 
limit on the nonthermal flux in the 0.2-10 keV range, L"^"^^^ < 0.08 x 10'^*' erg s~^, as found 
for PSR J2124-3358 in the two-temperature PC interpretation (§3.2). On the other hand, 
one can assume that the above-mentioned PC model is in fact the PC "core" and estimate 
temperature and size of the PC "rim" adding second thermal component. We found that to 
obtain a "rim" radius compatible with the estimates found for PSRs J0437-4715 and J2124- 
3358, a larger hydrogen density or smaller "core" radius is required. For example, if the PC 
"rim" temperature and radius are T™ = 0.4 MK and -R™ =1.1 km, then the other three 
parameters can be nH,20 = 3, Tp™'^ = 2.1 MK and R^^^^ = 0.07 km, or nH,20 = 6, Tp°^'^ = 1.7 
MK and R'^J'^ = 0.15 km. We note that although there exist many such combinations of 
model parameters, the main contribution to the thermal luminosity, about 75%, comes from 
the softer ("core") component as determined in the one-temperature PC fit. 
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4.3. Timing analysis 

The pulsar radio ephemeris parameters from the ATNF catalog, /o = 193.71568669103 
Hz and / = —6.953 x 10'^^ s~^ (at MJD 51018.0), were used for searching pulsations in 
the X-ray flux detected from PSR J1024-0719 with the EPIC-pn instrument. An optimal 
column interval in the 1-D distribution of EPIC-pn counts (Fig. 10) was determined with 
the aid of the Z\ test (other harmonics with n> \ were found to be less significant). It gave 
Z\ — 36.1 for 854 source- plus-background counts extracted from the CCD columns # 38-39 
in the 0.3-2 keV range. Varying spin frequency in vicinity of /o resulted in a maximal value 
-^imax — 39.1 at a frequency differing by bf = 2.3 /iHz from /q. The relative deviation, 
Sf/f = 1.2 X 10~^, is still consistent with the EPIC-pn timing accuracy. This value Zf^^^^ 
corresponds to a signal detection at a 6a significance level (probability to obtain this value 
by chance in one trial is 3 x 10^^). 

Figure 12 presents the light curve of PSR J1024-0719 extracted in the 0.3-2 keV range. 
The pulse profile reveals a single broad peak per period, in agreement with the result of 
the Z^^ test. To determine the intrinsic source pulsed fraction, we first estimated how many 
photons of those 854 extracted belong to the pulsar. To do this, 1-D count distribution were 
constructed from EPIC-MOS photons detected in the EPIC-pn filed of view divided in two 
parts, with and without the pulsar emission (see caption to Fig. 10), similar to what was 
done for PSR J0437-4715. It gave an estimate of 34% on the fraction of pulsar photons in the 
total number of EPIC-MOS counts of E ^ 0.3-2 keV extracted in the CCD columns # 38- 
39. Assuming that the same fraction stays for the EPIC-pn photons^^ gives the intrinsic 
pulsed fraction of about 52%, with a rather large error. 

5. PSR J0034-0534 

It is the third fastest pulsar known, with a spin period P ~ 1.9 ms. This pulsar is in a 
circular binary system with a low-mass companion and 1.6 d orbital period. The standard 
estimates give its characteristic age r ~ 6 Gyr, rotational energy loss E ^ 3.0 x 10^"^ erg s~^, 
and a distance to the pulsar d ^ 0.54 kpc. In earlier years PSR J0034-0534 was observed 
in X-rays only with BeppoSAX (in July 1999 for about 100 ks), but the pulsar was not 
detected. 



^^This assumption may not be very accurate because of, e.g., different internal background in the EPIC- 
MOS and EPIC-pn instruments. 
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XMM-Newton was pointed at PSR J0034-0534 on 2002 June 16 (orbit 463) for 35.0 
and 32.1 ks effective exposures with the EPIC-MOS and EPIC-pn detectors, respectively, 
with medium filters and the same observational modes as in the observations of the other 
three pulsar described §§ 2-4. A 20' x 20' combined EPIC-MOSl and M0S2 image is shown 
in Figure 13. In this image an X-ray source is evident at a position which differs by only 
0'.'2 from the radio position of PSR J2124-3358. Unfortunately, in the EPIC-pn observation 
with the satellite position angle of 65.5° this X-ray source was completely buried under 
emission from the much brighter (by a factor of 10) object indicated with the label "S" in 
the image. Therefore, neither spectral nor temporal information on this source could be 
obtained from the EPIC-pn data (timing analysis of EPIC-pn counts extracted at the source 
position revealed no pulsations). The count rates of the source at the pulsar position are 
0.9 ± 0.5 and 1.3 ± 0.5 counts ks^^ in the EPIC-MOSl and M0S2 detectors, respectively. 
Despite the low (< 3o") formal significance of the X-ray detection, the proximity of this 
source to the radio position of PSR J0034-0534 makes it a likely X-ray counterpart of the 
pulsar. About 70 counts detected from this source with the two EPIC-MOS detectors at 
E <3 keV can be equally well fitted with a single PL spectrum of F ~ 2.5 or a single BB 
model of ~ 2.2 MK and ~ 0.05 km (for d = 0.54 kpc). Both these models yields an 
X-ray luminosity Lx ~ 0.4 x 10^° erg s~^ 



6. Discussion and conclusion 

The new XMM-Newton observations have enabled us to perform the spectral and timing 
analyses of the X-ray data collected from three pulsars, J0437-4715, J2124-3358, and J1024- 
0719. In addition, a possible X-ray counterpart of PSR J0034-0534 has been found, with 
an estimate on its X-ray flux. For the first three objects the XMM-Newton data have 
revealed strong evidences of presence of thermal components in their X-ray emission. These 
components are interpreted as radiation emitted from heated PCs around magnetic poles on 
the NS surface. 

In the case of PSR J0437-4715, the new data confirmed the main result obtained from 
the previous observations that the pulsar spectrum is best fitted with a two-temperature 
PCs model with nonmagnetic hydrogen atmosphere plus a nonthermal (PL) component 
which contributes mainly at energies E >3 keV. In this interpretation, the PC temperature 
changes from about 1.5 MK in the PC "core" of a ~ 0.4 km radius down to ~ 0.5 MK in 
the PC "rim" of a ~ 2.5 km radius. The XMM-Newton data allowed us to investigate the 
energy-resolved pulse profiles in a broad range of = 0.3-6 keV, as well as perform phase- 
resolved spectroscopy. This analysis showed that the pulse shape is different at lower and 
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higher energies (below and above ~ 2 keV), with pulsed fraction /p increasing by about 10% 
in the 0.3-2 keV range where the thermal emission dominates. This increase of /p can be 
explained by the anisotropic properties of radiation emitted from a NS hydrogen atmosphere 
(see Zavlin et al. 1996, and Zavlin & Pavlov 1998 for more details). The pulsed fraction in the 
light curve extracted at energies above 2 keV is rather large, /p ^ 60%. In addition, the peak 
in the pulse profile obtained at higher energies is narrower than the pulse in the low-energy 
light curve and shifted by about 0.1 in phase. These properties provide strong evidence that 
two different emission mechanisms generate the X-ray emission of PSR J0437-4715. 

The XMM-Ncwton spectral data on PSR J2124-3358 can be most plausibly interpreted 
as either a combination of a one-temperature PC model plus a nonthermal PL spectrum or a 
pure thermal radiation from PCs with a nonuniform temperature decreasing from ~ 2 MK to 
~ 0.5 MK in the "core" and "rim" areas of the PCs, similar to that derived for PSR J0437- 
4715. Comparing the PC radii inferred in these two approaches — about 0.3 km with the 
PC-plus-PL combination and 1-2 km in the two-temperature PC model — with the canonical 
estimate of pulsar models, R^^ = R{2nR/cPy/^ ~ 2.1 km for PSR J2124-3358 (assuming 
the NS radius R = 10 km), suggests that the pure thermal interpretation is preferable. In 
any case, the thermal luminosities derived in these model fits differ insignificantly, by only 
a factor of ~ 2. In addition, the new data allowed us to firmly detect pulsations of the 
X-ray flux of PSR J2124-3358 at the radio pulsar period, although the shape of the obtained 
pulse profile cannot be interpreted unambiguously. It may be described as a profile with two 
peaks of different heights separated by about 0.35 (or 0.65) in phase. If so, it would suggest 
that two different components are detected in the pulsar emission. On the other hand, it 
could be a broad single peak produced by thermal PC emission with a shape distorted by 
the relativistic effects in fast rotating NSs. 

For PSR J2124-3358, the XMM-Newton observation provided the first firm detection 
of the pulsar X-ray emission. Although the obtained statistics is much poorer than that 
dehvered in the observations of PSR J2124-3358 and, especially PSR J0437-4715, it has 
shown that the pulsar's spectrum is best described as a pure thermal radiation emitted 
from PCs with either uniform or nonuniform ( "core" + "rim" ) temperature. The derived 
total thermal luminosity is almost independent on assumption about the PC temperature 
distribution. In addition to the spectral analysis, the new data revealed pulsations of the 
pulsar X-ray fiux at a highly significant level. Most likely, the pulsar light curves shows a 
single broad pulse per rotational period, although much better statistics is required to draw 
more definitive conclusions on the pulse shape. 

Regarding the XMM-Newton observation of PSR J0034-0534, only a possible X-ray 
counterpart of the pulsar could be suggested based on the very close proximity of this source 
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to the pulsar position. A longer observation, with a proper satellite position angle, is needed 
to perform a meaningful investigation of X-ray radiation of this pulsar. 

With a sample of millisecond pulsars showing thermal PC emission available (see Ta- 
ble 2), it is natural to compare observational properties of the PC radiation, first of all, the 
luminosities, with those predicted by pulsar models. A most recent and detailed analysis 
of PC heating is presented in Harding & Muslimov (2001, 2002) who modified the space 
charge limited flow model by Arons (1981) to account for the relativistic effects and elec- 
tric field screening above the NS surface. In these models the PC heating is produced by 
inverse Compton scattering (ICS) of thermal X-rays from the NS surface by primary elec- 
trons accelerated in the pulsar magnetosphere, as well as by positrons generated through 
curvature radiation (CR) and returning to the NS surface from the upper pair formation 
front. These authors showed that the PC heating from the CR pair fronts dominates for 
most pulsars except those with very short spin periods and low magnetic fields (see Fig. 1 in 
Harding & Muslimov 2002) because the latter objects cannot produced CR pairs. Therefore, 
the ICS is the only mechanism to heat PCs of millisecond pulsars. Harding & Muslimov 
(2002) calculated the "PC efficiency" of one PC, r]'^^ = L^'^i/E, for typical parameters of 
millisecond pulsars (age r = 0.1-10 Gyr, spin period P = 2-5 ms). This model predicts 1]^'^ 
ranging between about 5 x 10^^ and 5 x 10^^ for PCs with temperatures Tp^ = 1-3 MK 
(see Fig. 8 in Harding & Mushmov 2002). Comparing these estimates with those derived 
from the observations (Table 2) shows that the maximum predicted value of T^^f^x is lower 
by a factor of 10 than the results obtained for PSRs J0437-4715 and J0030+0451, and a 
factor of 6 for PSR J2124-3358. This difference is larger for PSR J1012+5307, but we note 
that the 1]^'^ value given in Table 2 for this pulsar was obtained under assumption that all 
X-ray flux detected from PSR J1012+5307 is of a thermal origin. For PSR J1024-0719, the 
value of ^7max "^^11 agrees with the observational result. On the other hand, the predicted PC 
luminosities, L^^ ~ lO^^-lO''^^ erg s"^ for Tpc = 1 MK and up to lO''^^ erg s^^ for Tpc = 3 MK 
(see Fig. 9 in Harding & Muslimov 2002), are in much better agreement with those estimated 
for the objects listed in Table 2. To conclude, the PC heating from ICS pairs is a promising 
model for interpreting properties of thermal PC emission observed from millisecond pulsars. 

Finally, as mentioned in § 1, there is no clear understanding yet at which conditions the 
thermal PC component is expected to dominate over the nonthermal one in X-ray radiation of 
millisecond pulsars. Saito et al. (1997) first suggested that high luminosity of the nonthermal 
emission may be associated with large values of the magnetic field at the pulsar light cylinder, 

— Bs^ri[R/ Ric]^ {Ric — cP/ 27r). Indeed, the four pulsars emitting pure nonthermal 
emission, PSRs B1937+21, B1957+20, B1821-24, and J0218+4232 (see §1 and Table 1), 
possess magnetic fields, Bic ~ (0.3-1) x 10® G, close to that of the Crab pulsar (~ 1 x 10® 
G) and exceeding those in the millisecond pulsars with thermal emission, B\c ~ (2-3) x 10^ 
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G, at least by an order of magnitude. This provides a strong indication that the nonthermal 
emission of millisecond pulsars is generated in emission zone(s) close to the light cylinder 
(as predicted by the outer-gap pulsar models — see, e.g., Cheng et al. 1986) and Bi^ is one 
of the main parameters to govern the magnetospheric activity. 

The author thanks George Pavlov for helpful and stimulating discussions. This work is 
supported by a National Research Council Research Associateship Award at NASA MSFC. 
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Table 1. Properties of nine millisecond pulsars^ 



PSR 


P 


d 


r 


logE 


logLx 


Refs. 




(ms) 


(kpc) 


(Gyr) 


(crgs-i) 


(crgs-i) 




B1937+21 


1.56 


3.57 


0.24 


36.04 


33.15 


1,2 


B1957+20 


1.61 


2.49 


2.24 


35.20 


31.81 


3 


J0218+4232 


2.32 


2.67 


0.48 


35.38 


32.54 


4 


B1821-24 


3.05 


3.09 


0.03 


36.34 


32.71 


5 


J0751+1807 


3.48 


1.15 


7.08 


33.86 


30.84 


6 


J0030+0451 


4.87 


0.32 


7.71 


33.53 


30.40 


7 


J2124- 3358 


4.93 


0.27 


3.80 


33.83 


30.23 


8 


J1012+5307 


5.26 


O.'ll 


4.86 


33.67 


.30.. 38 


6 


.10 137 1715 


5.7() 


D.ll 


l.S!) 


:->:-!. 5 S 


:-!0. l(i 


9 



^ Based on results obtained in previous works (see references below). 

Note. — Second through fourth columns give standcird pulsar parameters: 
spin periods P, distances d, characteristic ages r, and rotation energy losses E. 
All distances arc estimated from the pulsar dispersion measures and the model 
of Galactic distribution of free electrons (Cordcs & Lazio 2003), except for 
PSR J0437— 4715 with the distance determined from the pulsar's parallax (van 
Straten et al. 2001). The fifth column provides pulsars' X-ray luminosities 
in the 0.2—10 keV range as derived in the corresponding references (sixth 
column): 1 - Takahashi ct al. (2001), 2 - Nicastro ct al. (2004), 3 - Stappers 
et al. (2003), 4 - Webb et al. (2004a), 5 - Becker et al. (2004), 6 - Webb et 
al. (2004b), 7 - Becker & Aschenbach (2002), 8 - Sakurai et al. (2001) 9 - 
Zavlin et al. (2002) 
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Table 2. Thermal luminosities of millisecond pulsars 



PSR 


logB 
(ergs-i) 


(ergs 1) 


log JjP"^ 


J0437-4715 


33.58 


30.23 


-3.35 


J2124-3358 


33.83 


29.96 


-3.87 


J1024-0719 


33.72 


29.27 


-4.45 


J1012+5307'' 


33.67 


30.52 


-3.15 


J0030+0451 


33.53 


30.18 


-3.35 



Note. — The thermal (bolometric) luminosi- 
ties of one PC are determined as i^^i = ^bol/^) 
with Lbol estimates obtained in this work for the 
first three pulsars. For PSRs J0030-I-0451 and 
J1012-)-5307 the PC luminosities arc calculated as 
-'^bol — ^i^^^bol/^' "here Lj^i arc luminosities de- 
rived from BB spectral fits by Becker & Ascheubach 
(2002) and Webb et al. (2004b), respectively, and 
ffr = [1 - 2GM/c2il]-V2 (= 0.77 for M = 1.4 Mq 
and R = 10 km). The "PC efficiency" is defined as 
V"^ = LZi/E. 

^ Assuming that all flux detected for this pulsar 
is of a thermal origin. 
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Fig. 1.— Combined EPIC-MOSl and M0S2 image in the 0.3-5 keV range, with PSR J0437- 
4715 marked with the label "PSR". Long-dashed rectangle depicts the field of view of the 
EPIC-pn CCD used in the XMM-Newton observation of the pulsar. Arrow shows the direc- 
tion of the CCD read-out perpendicular to which EPIC-pn counts were collapsed in a 1-D 
distribution (see Fig. 2). Short-dashed line separates the areas within the EPIC-pn field of 
view with negligible contributions from the pulsar and AGN emission (above and below the 
line, respectively). 
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Fig. 2. — Upper panel: 1-D distribution of EPIC-pn counts in the 0.3-5 keV range. Peaks 
labeled with "PSR" and "AGN" are the images of PSR J0437-4715 and the neighboring 
AGN (see Fig. 1 and text). One CCD column has a width of 4'.'2. Lower panel: 1-D image 
(solid line) constructed from the combined EPIC-MOSl and M0S2 data confined within the 
EPIC-pn filed of view drawn in Fig. 1. Dotted (labeled with "PSR" ) and dashed (labeled with 
"AGN") histograms show 1-D EPIC-MOS images computed for the lower and upper areas 
of the EPIC-pn field of view (below and above the short-dashed line in Fig. 1, respectively). 
Thin vertical lines in the both panels indicate the CCD column below which contamination 
of the pulsar emission by that from the AGN is negligible. 
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Fig. 3.— EPIC-pn and two EPIC-MOS phase-integrated spectra of PSR J0437-4715 fit- 
ted witfi a best two-temperature polar cap (PC) model ( "core" -|- "rim" ) plus a power- low 
(PL) component, residuals of the fit, and model fluxes {upper, middle, and lower panels, 
respectively) . 



-24- 




Fig. 4. — Energy-resolved pulsed profiles of PSR J0437-4715 extracted from the EPIC-pn 
data, with estimated values of the intrinsic source pulsed fraction /p. 
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Fig. 5.— Pulsed profiles of PSR J0437-4715 obtained from the XMM-Newton (0.3-6 keV), 
Chandra (0.1-10 kcV) and ROSAT (0.1-2.4 kcV) data, with estimated values of the intrinsic 
source pulsed fraction /p. Zero phases are arbitrary. 
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Fig. 6. — Phase dependences of the polar cap (PC) areas and nonthermal (PL) flux in the 
two-temperature PC ( "core" + "rim" ) plus PL model. 
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Fig. 7.— Combined EPIC-MOSl and M0S2 image in the 0.3-3 keV range, with PSR J2124- 
3358 marked with the label "PSR" (see caption to Fig. 1 for other details). One- dimensional 
distribution of EPIC-pn counts is shown in the insert. 
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Fig. 8.— EPIC-pn and two EPIC-MOS phase-integrated spectra of PSR J2124-3358 fitted 
with a two-component (polar cap [PC] plus a power-low [PL]) model, residuals of the fit, 
and model fiuxes {upper, middle, and lower panels, respectively). 




Fig. 9.— Pulse profile of PSR J2124-3358 extracted from the EPIC-pn data in the 0.3-2 
keV range, with estimated value of the intrinsic source pulsed fraction /p. 
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Fig. 10.— Combined EPIC-MOSl and M0S2 image in the 0.3-5 keV range, with PSR J1024- 
0719 indicated with dotted arrow (see caption to Fig. 1 for other details). Emission of the 
pulsar dominates in the area between two short-dashed lines within the EPIC-pn field of 
view. One-dimensional distribution of EPIC-pn counts is shown in the insert. 
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Fig. 11.— Two EPIC-MOS phase-integrated spectra of PSR J1024-0719 fitted with a polar 
cap model, and residuals of the fit. 
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Fig. 12.— Pulse profile of PSR J1024-0719 extracted from the EPIC-pn data in the 0.3-2 
keV range, with estimated value of the intrinsic source pulsed fraction /p. 
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Fig. 13.— Combined EPIC-MOSl and M0S2 image in the 0.3-3 keV range, with PSR J0034- 
0534 indicated with dotted arrow (see caption to Fig. 1 for other details). The brightest 
source in the filed of view is marked with the label "S" . 



